A Nonadditive Quantum Code 

Eric M. Rains, R. H. Hardin, Peter W. Shor, and N. J. A. Sloane 
AT&T Research, 600 Mountain Ave., Murray Hill, NJ 07974, USA 
(March 3, 1997) 



On 
OS 



(N ' 

> : 

(N . 
O ■ 

o ■ 
cn : 
o . 
r» ■ 

9 L,: 
^— > ■ 

c ■ 

3 : 
cr 



X 



Up to now every good quantum error-correcting code dis- 
covered has had the structure of an eigenspace of an Abelian 
group generated by tensor products of Pauli matrices; such 
codes are known as stabilizer or additive codes. In this letter 
we present the first example of a code that is better than any 
code of this type. It encodes six states in five qubits and can 
correct the erasure of any single qubit. 



Most present models for quantum computers require 
quantum error-correcting codes for their operations. 
Many recent papers have been devoted to the construc- 
tion of these codes |^-||] . Up to now all good codes known 
have fallen into the class of what have been called "sta- 
bilizer" or "additive" codes @,§. These are not the most 
general quantum codes, however, and it was an open 
question as to whether more general codes could be more 
powerful. In the present letter we present a code which 
is not a stabilizer code, and which is strictly better than 
any code of this type. 

Stabilizer codes are defined as follows. Let the ex- 
traspecial group E be generated by all tensor prod- 
ucts Mi <g> M2 ® . . . ® M n , where each Mi is one of 
{I,a x ,a y ,a z }, where 



I = 



-0) 



a~ x , o~ y , <j z being the Pauli matrices. Then a stabilizer 
code is a joint eigenspace of an Abelian subgroup of E. 

The new code. Since a quantum code is a subspace of 
the Hilbert space C 2 , we may define it by giving the or- 
thogonal projection matrix onto this subspace. Since the 
Pauli matrices together with I form a basis for the space 
of 2 x 2 Hcrmitian matrices, we may write the projection 
as a linear combination of elements of E. In this form, 
the new code Q is defined by the projection matrix 

P = 1/16 [3 I®I®I®I®I 

+ (I ® a z ® a y ® a y ® a z ) cyc 
+ (I<E) a x ®C7 Z <g>er z ®a x ) cyc 
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where the subscript "eye" indicates that all five cyclic 
shifts occur. It is straightforward to verify that P 2 = P 



and Trace(P) = 6, so that the eigenvalues of P are 1 
(6 times) and (26 times). Thus P projects onto a six- 
dimensional space. 

It remains to verify that the code has minimum dis- 
tance 2, that is, that it can correct one erasure. For a 
code to correct one erasure it suffices for it to be orthog- 
onal to its image under any single qubit error. Since Q is 
cyclic, it is enough to consider an error in the first qubit. 
There are thus three cases to check: conjugation of P 
by r) = M®I®I<g>I®I®I where M is one of a x , 
o~y , <j z , obtaining P' = rjPr). We verify in each case that 
PP' = 0. Thus Q is a ((5,6,2)) code, as claimed. 

Further properties of the code. Since the code is 
cyclic, we label the qubits by 0, 1, ... ,4 (mod 5). The 
code also has the symmetry defined by 



k — > 2k (mod 5), a x — > a y 



(1) 



There are further symmetries that do not permute the 
qubits, for example conjugation by the element 

<y z ® cr x (8) a y ® a y (X> o~ x (2) 

of E, or any of its cyclic shifts. These five elements 
generate a subgroup of E of size 32, which we call H . 
Moreover, since this group is Abelian, it has 32 one- 
dimensional eigenspaces, and Q is the span of six of them. 
(Each of these eigenspaces is a ((5,1,3)) stabilizer quan- 
tum code.) This gives an explicit basis for our code, 
namely: 

looooo) - (|oooii»cyc + (|ooioi)) cyc - (|omi)) cyc (3) 

together with all five cyclic shifts of 

|00001) - |00010) - |00100) - |01000) - |10000) 
+ |00111) - |01110) - |11100) + |11001) + |10011) 
- |01011) + 10110) - |01101) + |11010) - |10101) 
-|11111) (4) 

By examining how single bit errors act on these 32 
eigenspaces, one can again see that the image of Q under 
such an error is orthogonal to Q, by noting that the six 
new eigenspaces are disjoint from the original six. 

The symmetries discussed above generate a group of 
order 640. There is an additional symmetry, which can 
be described as follows: First, permute the columns as 
k — > k 3 , that is, exchange the qubits 2 and 3. Next, 
for each qubit, negate one of the Pauli matrices and ex- 
change the other two, where the Pauli matrices negated 



1 



are er 2 , a y , a x , cr x , <7 y , respectively. This increases the 
size of the symmetry group to 3840. The group acts as 
the permutation group 5s on the qubits. This is the full 
group of symmetries of the code, that is, the full sub- 
group of the semidirect product of S5 and PSU(2) 5 that 
preserves the code [jl0| . 

In principal one can construct other codes in a similar 
manner, e.g. as the spans of translates of self-dual stabi- 
lizer codes. Let Co be a self-dual additive code of length 
n with associated stabilizer quantum code Qo and let C 
be the union of K cosets of Cq. If C has minimum dis- 
tance d, then there exists an ((n, K, d)) quantum code. 
This code is obtained as the span of the translates of Qo 
under the operators in C. For our code, for instance, Co 
is the [[5,0,3]] additive code H and C is Co + (cCo) cy c! 
where 

<j = Gt.®(j 7 ®I®I®I. 



Thus, for example, we have 

P = P + (aPoOcyc, 



(5) 



where Pq is \x)(x\ and \x) is given by Eq. (3). 

We remark that the only stabilizer code that contains 
Q is the trivial code consisting of the entire Hilbert space. 

How the code was discovered. Investigation of the lin- 
ear programming bound for quantum codes (cf. |||ll||l2| ) 
revealed that for a length 5 code of minimum distance 2 
the dimension K = 6 is extremal. In fact, the inequalities 
cannot even be satisfied for any real number K > 6. This 
suggested a ((5,6,2)) code as a natural place to look. Fur- 
ther, the weight enumerator of any such code is uniquely 
determined by linear programming; with the normaliza- 
tion used in [O], it is 



A(u, v) = 36u 5 + 60uv 4 + 96v 5 



(6) 



The best stabilizer code of this length and distance has 
dimension K = 4 || Table III]. 

Rather than look for a projection matrix with this 
weight enumerator, we decided to look for a matrix with 
this weight enumerator that happened to be a projec- 
tion. We describe a near projection matrix M as a linear 
combination of elements of E. Beginning with a random 
matrix M, we first scaled the contribution from the terms 
of each weight so as to force the desired weight enumer- 
ator. We then replaced M by M' = 2M 2 - M 4 . (This 
moves small eigenvalues of M closer to and large eigen- 
values closer to 1.) We then iterated these two steps. On 
our first attempt, the algorithm converged to a (complex) 
projection matrix. Since the algorithm preserves symme- 
try, by choosing our initial guess to be real, we obtained 
a real code equivalent to the code described above. We 
omit the details of transforming the computer output into 
the above form. Note that we did not choose our initial 
guess to be cyclic, but nonetheless obtained a code equiv- 
alent to one with cyclic symmetry. 
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